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Abstract 

n— t ' Quantum magnetism is one of the most active areas of research in con- 

pH , densed matter physics. There is significant research interest specially in 

Q ' low-dimensional quantum spin systems. Such systems have a large num- 

O ', ber of experimental realizations and exhibit a variety of phenomena the 

origin of which can be attributed to quantum effects and low dimensions. 

In this review, an overview of some aspects of quantum magnetism in 

►^ [ low dimensions is given. The emphasis is on key concepts, theorems and 

Qs . rigorous results as well as models of spin chains, ladders and frustated 

O^ ' magnetic systems. 

m ■ 

O ■ 1 Introduction 

Quantum magnets are spin systems in which the spins interact via the well- 
_. , known exchange interaction. The interaction is purely quantum mechanical in 

e nature and the form of the interaction was derived simultaneously by Heisenberg 

I and Dirac in 1926 |Q]. The most well-known model of interacting spins in an 

n^ ■ insulating soHd is the Heisenberg model with the Hamiltonian 

O. h^J2JvS.-S, (1) 



Si is the spin operator located at the lattice site i and Jy denotes the strength 



S^ 



can have a magnitude 1/2, 1, 3/2, 2, 



of the exchange interaction. The spin 
5^ , ...etc. The lattice, at the sites of which the spins are located, is d-dimensional. 

Examples are a Hnear chain (d = 1), the square lattice (d = 2 ) and the cubic 
lattice (d =3 ) . Ladders have structures interpolating between the chain and the 
square lattice. An n-chain ladder consists of n chains (n = 2, 3, 4, ...etc.) coupled 
by rungs. Real magnetic solids are three-dimensional (3d) but can be effectively 
considered as low-dimensional systems if the exchange interactions have different 



strengths in different directions. To give an example, a magnetic solid may 
consist of chains of spins. The solid may be considered as a linear chain (d=l) 
compound if the intra-chain exchange interactions are much stronger than the 
inter-chain ones. In a planar (d=2) magnetic system, the dominant exchange 
interactions are intra-planar. Several examples of low-dimensional magnetic 
systems are given in |p|. 

The strength of the exchange interaction Jij in Eq.(l) falls down rapidly 
as the distance between interacting spins increases. For many solids, the sites 
i and j are nearest-neighbours (n.n.s) on the lattice and Jij 's have the same 
magnitude J for all the n.n. interactions. The Hamiltonian in (1) then becomes 

H = JJ2S:.S', (2) 

There are, however, examples of magnetic systems in which the strengths of the 
exchange interactions between successive pairs of spins are not the same. Also, 
the interaction Hamiltonian (1) may include n.n. as well as further-neighbour 
interactions. The well-known Majumdar-Ghosh chain |g| is described by the 
Hamiltonian 

N J N 

Hmg = J 2_^ Si.S i+i + — 2_^ Si . S i+2 (3) 

i=l i=l 

and includes n.n. as well as next-nearest-neighbour (n.n.n.) interactions. The 
Haldane-Shastry model U] has a Hamiltonian of the form 

and includes long-ranged interactions. Real materials are characterised by var- 
ious types of anisotropy. The fully anisotropic n.n. Heisenberg Hamiltonian in 
Id is given by 

N 

HxYZ = ^JiJxSiS"^^^ + JySfSf^^ + J:,SiS-^^] (5) 

1=1 

The special cases of this Hamiltonian are: the Ising {Jx = Jy = 0) , the XY 
{Jz = 0), the XXX or isotropic Heisenberg {Jx = Jy — Jz ) and the XXZ or 
anisotropic Heisenberg {Jx = Jy ^ Jz ) models. There is a huge literature 
on these models some of which are summarised in Refs. ||^, ||, ||, ^, Q- Other 
anisotropy terms may be present in the full spin Hamiltonian besides the basic 
exchange interaction terms. 

Consider the isotropic Heisenberg Hamiltonian in (2) where (zj) denotes 
a n.n. pair of spins. The sign of the exchange interaction determines the 
favourable alignment of the n.n. spins. J > 0( J < 0) corresponds to anti- 
ferromagnetic (ferromagnetic) exchange interaction. To see how exchange in- 
teraction leads to magnetic order, treat the spins as classical vectors. Each n.n. 



spin pair has an interaction energy JS^cos9 where 9 is the angle between n.n. 
spin orientations. When J is < 0, the lowest energy is achieved when 6* = 0, 
i.e., the interacting spins are parallel. The ferromagnetic (FM) ground state 
has all the spins parallel and the ground state energy Eg — ~J^^^- where z 
is the coordination number of the lattice. When J is > 0, the lowest energy is 
achieved for 9 = tt , i.e., the n.n. spins are antiparallel. The antiferromagnetic 
(AFM) ground state is the Neel state in which n.n. spins are antiparallel to 
each other. The ground state energy Eg = — ^^^^ . 

Magnetism , however, is a purely quantum phenomenon and the Hamiltonian 
(2) is to be treated quantum mechanically rather than classically. For simplic- 
ity, consider a chain of spins of magnitude i. Periodic boundary condition is 
assumed, i.e., S n+i = S\ . The Hamiltonian (2) can be written as 

^ 1 

H = JJ2[S:SU, + -{S+Sr^, + 5r5+ i)] (6) 

i=l 

where 

St=Sf±^Sf (7) 

are the raising and lowering operators. It is easy to check that in the case of a 
FM, the classical ground state is still the quantum mechanical ground state with 
the same ground state energy. However, the classical AFM ground state (the 
Neel state) is not the quantum mechanical ground state. The determination of 
the exact AFM ground state is a tough many body problem and the solution 
can be obtained with the help of the Bethe Ansatz technique (Section 2). 

For a spin-1/2 system, the number of eigenstates of the Hamiltonian is 2^ 
where N is the number of spins. In a real solid N is ~ 10^"^ and exact determi- 
nation of all the eigenvalues and the eigenfunctions of the system is impossible. 
There are some classes of AFM spin models for which the ground state and in 
a few cases the low-lying excitation spectrum are known exactly (Section 2) . In 
the majority of cases, however, the ground state and the low-lying excited states 
are determined in an approximate manner. Knowledge of the low-lying excita- 
tion spectrum enables one to determine the low-temperature thermodynamics 
and the response to weak external fields. The usual thermodynamic quantities 
of a magnetic system are magnetisation, specific heat and susceptibility. Ex- 
change interaction can give rise to magnetic order below a critical temperature. 
However, for some spin systems, the ground state is disordered, i.e., there is no 
magnetic order even at T = 0. Long range order (LRO) of the Neel-type exists 
in the magnetic system if 

UmR^oo ( 5^(0). S^(i?)) ^ (8) 

where R denotes the spatial location of the spin. At T = 0, the expecta- 
tion value is in the ground state and at T 7^ , the expectation value is the 
usual thermodynamic average. The dynamical properties of a magnetic system 



are governed by the time-dependent pair correlation functions or their Fourier 
transforms. Quantities of experimental interest include the dynamical corre- 
lation functions in neutron scattering experiments, the NMR spin-lattice re- 
laxation rate, various relaxation functions and associated Hneshapes as well as 
the dynamical response of the magnetic system to various spectroscopic probes 
[|j . Knowledge of the ground and low- lying states and the corresponding energy 
eigenvalues is essential to determine the thermodynamic and dynamic properties 
of a magnetic system. 

The discovery of high-temperature cuprate superconductors in 1987 has 
given a tremendous boost to research activity in magnetism. The dominant 
electronic and magnetic properties of the cuprate systems are associated with 
the copper-oxide {CUO2) planes. The Cv?^ ions carry spin-i and the spins 
interact via the Heisenberg AFM exchange interaction Hamiltonian. This fact 
has given rise to a large number of studies on 2d antiferromagnets. The cuprates 
exhibit a variety of novel phenomena in their insulating, metallic and supercon- 
ducting phases some of which at least have links to quantum magnetism. The 
subject of magnetism has, as a result, expanded significantly in scope and con- 
tent. A rich interplay between theory and experiments has led to the discovery 
of materials exhibiting hitherto unknown phenomena, formulation of new theo- 
retical ideas, solution of old puzzles and opening up of new research possibilities. 
In this review, a brief overview of some of the important developments in quan- 
tum magnetism will be given. The focus is on quantum antiferromagnets and 
insulating solids. 

2 Theorems and rigorous results 

(i) Theorems : 

A. Lieb-Mattis theorem |y 

For general spin and for all dimensions and also for a bipartite lattice, the 
entire eigenvalue spectrum satisfies the inequality 

Eo{S)<E^{S + l) (9) 

where Eq(S) is the minimum energy corresponding to total spin S. The weak 
inequality becomes a strict inequality for a FM exchange coupHng between spins 
of the same sublattice. The theorem is valid for any range of exchange coupling 
and the proof does not require PBC. The ground state of the S — ^ Heisenberg 
AFM with an even number N of spins is a singlet according to the Lieb-Mattis 
theorem. 

B. Marshall's sign rule Q 

The rule specifies the structure of the ground state of a n.n. S = ^ Heisen- 
berg Hamiltonian defined on a bipartite lattice. The rule can be generaHsed to 
spin S, n.n.n. FM interaction but not to n.n.n. AFM interaction. A bipartite 
lattice is a lattice which can be divided into two sublattices A and B such that 
the n.n. spins of a spin belonging to the A sublattice are located in the B sub- 
lattice and vice versa. Examples of such lattices are the Hnear chain, the square 



and the cubic lattices. According to the sign rule, the ground state ip has the 
form 

\^) = Y,c^\^^) (10) 

A* 

where |/i) is an Ising basis state. The coefficient C^ has the form 

C, = {-r^a, (11) 

with a^ real and > and p^ is the number of up-spins in the A sublattice. 

C. Lieb, Schultz and Mattis (LSM) theorem [|ll|: 

A half-integer S spin chain described by a reasonably local Hamiltonian re- 
specting translational and rotational symmetry either has gapless excitation 
spectrum or has degenerate ground states, corresponding to spontaneously bro- 
ken translational symmetry. 

In the case of a gapless excitation spectrum, there is at least one momentum 
wave vector for which the excitation energy is zero. For a spectrum with gap, 
the lowest excitation is separated from the ground state by an energy gap A. 
The temperature dependence of thermodynamic quantities is determined by 
the nature of the excitation spectrum (with or without gap). The LSM theorem 
does not hold true for integer spin chains. For such chains, Haldane made a 
conjecture that the spin excitation spectrum is gapped [l^ . This conjecture 
has been verified both theoretically and experimentally |13|| . 

D. Oshikawa, Yamanaka and Affieck theorem ||lj| 

This theorem extends the LSM theorem to the case of an applied magnetic 
field. The content of the theorem is : translationally invariant spin chains in an 
applied field can have a gapped excitation spectrum, without breaking transla- 
tional symmetry, only when the magnetization per site m. {m — j^ J2i=i ^t ^ ^ 
is the total number of spins in the system ) obeys the relation 

S — m ~ integer (12) 

where S is the magnitude of the spin. The proof is an easy extension of that of 
the LSM theorem. The gapped phases correspond to magnetization plateaux in 
the m vs. H curve at the quantized values of m which satisfy (12). Whenever 
there is a gap in the spin excitation spectrum, it is obvious that the magnetiza- 
tion cannot change in changing external field. Fractional quantization can also 
occur, if accompanied by (expHcit or spontaneous) breaking of the translational 
symmetry. In this case, the plateau condition is given by 

n{S — m) = integer (13) 

where n is the period of the ground state. Hida ||lj has considered & S — \ 
HAFM chain with period 3 exchange coupHng. A plateau in the magnetization 
curve occurs at to = ^ (i of full magnetization ). In this case, n =3, S = ^ and 
m = ^ and the quantization condition in (13) is obeyed. Ref. |Q gives a review 
of magnetization plateaux in interacting spin systems. Magnetization plateaux 



have been observed in the magnetic compound NH4CUCI3 at m = j and | 
[ p^ . Possible extensions of the LSM theorem to higher dimensions have been 
suggested ||l^. The compound SrCu2(BO^)2 is the first AFM compound in 2d 
in which magnetization plateaux have been observed experimentally [^. Like 
the quantum Hall effect, the phenomenon of magnetization plateaux is another 
striking example of the quantization of a physically measurable quantity as a 
function of the magnetic field. 
E. Mermin- Wagner's theorem [g^ 

There cannot be any AFM LRO at finite T in dimensions d =1 and 2. The 
LRO can, however, exist in the ground state of spin models in d =2. LRO exists 
in the ground state of the 3d HAFM model for spin S > ^ |pi[. At finite T, the 
LRO persists upto a critical temperature Tc . For square p^ and hexagonal 
[ P3| lattices, LRO exists in the ground state for S* > 1 . The above results are 
based on rigorous proofs. No such proof exists as yet for S" = ^ , d =2 (this 
case is of interest because the Cu02 plane of the high- Tc cuprate systems is a 
5= i 2d AFM). 
(ii) Exact Results : 
A. the Bethe Ansatz ||i) 

The Bethe Ansatz (BA) was formulated by Bethe in 1931 and describes a 
wave function with a particular kind of structure. Bethe considered the spin— ^ 
Heisenberg linear chain in which only n.n. spins interact. In the case of the 
FM chain, the exact ground state is simple with all spins aligned in the same 
direction, say, pointing up. An excitation is created in the system by deviating 
a spin from its ground state arrangement, i.e., replacing an up-spin by a down- 
spin. Due to the exchange interaction, the deviated spin does not stay localised 
at a particular site but travels along the chain of spins. This excitation is the 
so-called spin wave or magnon. For the isotropic FM Heisenberg Hamiltonian, 
the exact one-magnon eigenstate is given by 

N 

V,= ^e^^"^„|m ) (14) 

m— 1 

where m denotes the site at which the down-spin is located and the summation 
over m runs from the first to the last site in the chain. The k's are the "momenta" 
which from periodic boundary conditions have N allowed values 

A:=^A,A = 0,l,2,...,iV-l (15) 

The excitation energy e^, measured w.r.t. the ground state energy and in units 
of J, is 

e/c = (1 — cosk) (16) 

In the case of r spin deviations (magnons) , the eigenfunction can be written as 

V'^''^ = Yl a(wim2,...,mr)5„^S'~^....S'„JTTT ) (17) 

7711 <m2<---'m.r 



The amplitudes are given by the BA 

a(mi,m2,...,m,) =^e'^.""^'"^ + ^'^.'<>"-"' (18) 

p 

where P stands for a permutation of the set {l,2,...,r} and Pj is the image of 
j under permutation. The sum is over all the r! permutations. Each term in 
(18) describes r plane waves. The scattering of a pair of waves introduces the 
phase shift ipji . The symmetric sum over permutations is in accordance with 
the bosonic nature of the waves, the spin waves, propagating along the chain. 
The energy of the state -0 is 

r 
i=l 

The k's are determined as before by applying the periodic boundary conditions 
which leads to the r equations 

Nk^ = 27rA, + Y^ (l)ij (20) 

3 

where Aj 's are r integers. One further imposes the condition that a spin at a 
particular site cannot be deviated more than once leading to the relations 

1 k k- 

2cot-(j)tj = cot-^ - cot-^ (21) 

Since cf)ij — ~(j)ji , Eqs. (21) are '^^^^ ' in number, i.e., there are as many 

r(r+l) 



2 



distinct 4> 's. Eqs. (20) are r in number. Together, the total number is 
equations in as many unknowns. Bethe thus established that the set of equations 
could be expected to have solutions. 

The momenta fci 's can be real or complex. In the first case, the spin waves 
or magnons scatter against each other giving rise to a continuum of scattering 
states. In the second case, the magnons form bound states, i.e., the reversed 
spins tend to be located at n.n. lattice positions. For r magnons, the r-magnon 
bound state energy is given by 

e=i(l-cosX) (22) 

r 

where K = X]i=i ^i is the total centre of mass momentum of r magnons. The re- 
sults can be generalised to the XXZ Heisenberg Hamiltonian. The multimagnon 
bound states were first detected in the quasi-one-dimensional magnetic system 
CoC^2-2i?20 at pumped helium temperatures and in high magnetic fields by 
far infrared spectroscopy [g5|. Later improvements ^^ made use of infrared 
HCN/DCN lasers, the high intensity of which made possible observation of 
even 14 magnon bound states. 

The exact ground state energy of the isotropic Heisenberg Hamiltonian 
(Eq.(2)) can be determined using the BA. The BA equations are the same 



Eg^ — JNln2 (24) 



as in the FM case but the sign of the exchange integral changes from —J to J 
(J > ). The total spin of the AFM ground state is S" = according to the 
Lieb-Mattis theorem. In the ground state, -^ spins are up and -y spins down 
(r = y)- The ground stae is non-degenerate and there is a unique choice of the 

Xi 's as 

Ai = 1,A2 = 3,A3 = 5, ....,A« =iV-l (23) 

The ground state is also spin disordered, i.e., has no AFM LRO. The exact 
ground state energy Eg is 

NJ 

T 

The low-lying excitation spectrum has been calculated by des Cloizeaux and 
Pearson (dCP) |2^ by making appropriate changes in the distribution of Xi 's 
in the ground state. The spectrum is given by 

TT 

e — — \sink\ , — tt < k < tt (25) 

for spin 1 states. The wave vector k is measured w.r.t. that of the ground 
state. A more rigorous calculation of the low-lying excitation spectrum has 
been given by Faddeev and Takhtajan [g8|. There are 5 = 1 as well as S* = 
states. We give a qualitative description of the excitation spectrum, for details 
Ref. ||2§] should be consulted. The energy of the low-lying excited states can be 
written as E{ki,k2) = e(fci) + £(^2) with e(fci) — ^sinki and total momentum 
k = ki + k2. At a fixed total momentum fc, one gets a continuum of scattering 
states. The lower boundary of the continuum is given by the dCP spectrum 
(one of the k'^s — ). The upper boundary is obtained for fci = fc2 = | and 



k 

sin— 

2 



(26) 



The energy-momentum relations suggest that the low-lying spectrum is actually 
a combination of two elementary excitations known as spinous. The energies 
and the momenta of the spinous just add up, showing that they do not interact. 
A spinon is a 5 = | object, so on combination they give rise to both 5" = 1 and 
S = states. In the Heisenberg model, the spinous are only noninteracting in 
the thermodynamic limit N ^(x. For an even number N of sites, the total spin 
is always an integer, so that the spins are always excited in pairs. The spinous 
can be visualised as kinks in the AFM order parameter. Due to the exchange 
interaction, the individual spinous get delocalised into plane wave states. In- 
elastic neutron scattering study of the linear chain S = ^ HAFM compound 
KCuFs has confirmed the existence of unbound spinon pair excitations p9| . 
The Haldane-Shastry model Q is another spin— ^ model in Id for which the 
ground state and low-lying excitation spectrum are known exactly. The ground 
state has the same functional form as the fractional quantum Hall ground state 
and is spin-disordered. The elementary excitations are spinous which are non- 
interacting even away from the thermodynamic limit, i.e., in finite systems. 



The individual spinons behave as semions, i.e., have statistical properties inter- 
mediate between fermions and bosons. In the case of integer spin chains, the 
spinons are bound and the excitation spectrum consists of spin- wave-like modes 
exhibiting the Haldane gap. The BA technique described in this Section is the 
one originally proposed by Bethe. There is an algebraic version of the B A which 
is in wide use and which gives the same final results as the earlier technique. For 
an introduction to the algebraic BA method, see the Refs. pO| , pl|. A tutorial 
review of the BA is given in Ref. [g2|- The BA was originally proposed for 
the Heisenberg model in magnetism. Later, the method was appHed to other 
interacting many body systems in Id such as the Fermi and Bose gas models 
in which particles on a line interact through delta function potentials ||33|, the 
Hubbard model in Id |Q, Id plasma which crystallizes as a Wigner solid ||35| , 
the Lai-Sutherland model l^^l which includes the Hubbard model and a dilute 
magnetic model as special cases, the Kondo model in Id |^, the single impurity 
Anderson model in Id ||3§], the supersymmetric t-J model (J = 2t) jsg] etc. In 
the case of quantum models, the BA method is applicable only to Id models. 
The BA method has also been applied to derive exact results for classical lattice 
statistical models in 2d. 
B. The Majumdar- Ghosh chain ||, 0| 

The Hamiltonian is given in Eq. (3) . The exact ground state of Hmg is doubly 
degenerate and the states are 

01 = [12][34]...[iV - liV], (/.2 = [23][45]...[iVl] (27) 

where [Im] denotes a singlet spin configuration for spins located at the sites 
I and m. Also, PBC is assumed. One finds that translational symmetry is 
broken in the ground state. The proof that 0i and 02 are the exact ground 
states can be obtained by the method of 'divide and conquer'. One can verify 
that 01 and 02 are exact eigenstates of Hmg by applying the spin identity 
Sn-(Si + S„,)[li^] = . Let El be the energy of 0i and 02 . Let Eg be the 
exact ground state energy. Then Eg < Ei . One divides the Hamiltonian H into 
sub-Hamiltonians , Hi 's, such that H — J^i Hi . Hi can be exactly diagonahsed 
and let EiQ be the ground state energy. Let ■0g be the exact ground state wave 
function. By variational theory. 

Eg = {^g \H\ ^Pg) ^ ^ (^g | iJ, | ^^ ) > ^ ^'^ (^S) 

i i 

One thus gets, 

J2 ^« <Eg<Ei (29) 

■i 

If one can show that J2i ^io — Ei , then Ei is the exact ground state energy. 
For the MG-chain, the sub-Hamiltonian Hi is 

Hi — —{ Si. S i+i + o i+i- o i+2 + <5 i+2- Si) (30) 



There are N such sub-Hamiltonians. One can easily verify that Eio — —^ and 
-El = —-XT ( "T ^^ ^^^ energy of a singlet and there are ^ VBs in (j)i and 
(^2)- From (29), one finds that the lower and upper bounds of Eg are equal and 
hence 0i and 02 are the exact ground states with energy Ei — —^^ ■ There 
is no LRO in the two-spin correlation function in the ground stae: 

K\^,J) = {S^S^)^\s,,-l6^,_^^^, (31) 

The four-spin correlation function has off-diagonal LRO. 



+ ^^\i-j\,i^\i~m\Aexp{iT:{— —)) (32) 

Let T be the translation operator for unit displacement. Then 

T(j)i = (j)2,T(f)2 = <j)i (33) 

The states 

0+ = -i=(01+ 02), r = -^(01-02) (34) 

correspond to momentum wave vectors A: = and k — n . The excitation 
spectrum is not exactly known. Shastry and Sutherland jio) have derived the 
excitation spectrum in the basis of 'defect' states. A defect state has the wave 
function 



ip{p,m) = ...[2p-3.2p-2]a2p-i[2p, 2p-|-l]...[2m-2,2TO-l]a2m[2m+l,2TO+2]... 

(35) 
where the defects {a2p-i and a2m) separate two ground states. The two defects 
are up-spins and the total spin of the state is 1. Similarly, the defect spins can 
be in a singlet spin configuration so that the total spin of the state is 0. Because 
of PEC, the defects occur in pairs. A variational state can be constructed 
by taking a linear combination of the defect states. The excitation spectrum 
consists of a continuum with a lower edge at J(| — 2 \cosk\). A bound state 
of the two defects can occur in a restricted region of momentum wave vectors. 
The MG chain has been studied for general values aJ of the n.n.n. interaction 
[ pj . The ground state is known exactly only at the MG point a = \ . The 
excitation spectrum is gapless for < a < acr{— 0.2411). Generalizations of 
the MG model to two dimensions exist |^, Q. The Shastry-Sutherland model 
[ ^ is defined on a square lattice and includes diagonal interactions as shown 
in Figure 1. The n.n. and diagonal exchange interactions are of strength Ji 
and J2 respectively. For jp- below a critical value ~ 0.7 , the exact ground state 
consists of singlets along the diagonals. At the critical point, the ground state 



10 



changes from the gapful disordered state to the AFM ordered gapless state. The 
compound SrCu2{BO^)2 is well-described by the Shastry-Sutherland model 
[ p^ . Bose and Mitra iQ have constructed a Ji — J2 — J3 — J4 — J5 spin-i 
model on the square lattice. Ji, J2, J3, J4 and J5 are the strengths of the n.n., 
diagonal, n.n.n., knight 's-move-distance-away and further-neighbour-diagonal 
exchange interactions (Figure 2). The four columnar dimer states (Figure 3) 
have been found to be the exact eigenstates of the spin Hamiltonian for the 
ratio of interaction strengths 

Ji : J2 : J3 : J4 : J5 = 1 : 1 : 2 : 2 ■ 4 ^^^^ 

It has not been possible as yet to prove that the four columnar dimer states are 
also the ground states. Using the method of 'divide and conquer', one can only 
prove that a single dimer state is the exact ground state with the dimer bonds 
of strength 7 J. The strengths of the other exchange interactions are as specified 
in (36). For a 4 x 4 lattice with PBC, one can trivially show that the four CD 
states are the exact ground states. 

C. The Affieck-Kennedy-Lieb-Tasaki (AKLT) model @ 
We have already discussed the LSM theorem, the proof of which fails for 
integer spin chains. Haldane | p2| , p^ in 1983 made the conjecture, based on 
a mapping of the HAFM Hamiltonian, in the long wavelength limit, onto the 
nonlinear a model, that integer-spin HAFM chains have a gap in the excitation 
spectrum. The conjecture has now been verified both theoretically and exper- 
imentally ||45J. In 1987, AKLT constructed a spin-1 model in Id for which the 
ground state could be determined exactly p4|. Consider a Id lattice, each site 
of which is occupied by a spin-1. Each such spin can be considered to be a 
symmetric combination of two spin-^ 's. Thus, one can write down 



V'-- = I— ),^^ = -i 

i>+- = ^(|+-> + !-+>, 5^ = 

^-+ = yj+- (37) 

where '+' ('— ') denotes an up (down) spin. 

AKLT constructed a valence bond solid (VBS) state in the following man- 
ner. In this state, each spin-^ component of a spin-1 forms a singlet (valence 
bond) with a spin-^ at a neighbouring site. Let e"^ {a,P = + or — ) be the 
antisymmetric tensor: 

£++ = f— = 0, e+- = -£-+ = 1 (38) 

A singlet spin configuration can be expressed as -y^^'^^ WP) , summation over 
repeated indices being implied. The VBS wave function (with PBC) can be 
written as 



11 



IV'VBs) is a linear superposition of all configurations in which each S^ — +1 is 
followed by a S*^ = — 1 with an arbitrary number of 5^ = spins in between 
and vice versa. If one leaves out the zero's, one gets a Neel-type of order. One 
can define a non-local string operator 

i-i 

a" = ~STexp{m V St)Sf, (a = x, y, z) (40) 



l=i+l 



and the order parameter 



0"t„„g = lim\i_^\^^ {af^ (41) 

The VBS state has no conventional LRO but is characterised by a non-zero 
value I of 0"(^j„g. After constructing the VBS state, AKLT determined the 
Hamiltonian for which the VBS state is the exact ground state. The Hamiltonian 

is 

Haklt = Y. ^2(5^* + ^^+l) (42) 

where P2 is the projection operator onto spin 2 for a pair of n.n. spins. The 
presence of a VB between each neighbouring pair impHes that the total spin 
of each pair cannot be 2 (after two of the S — i^ variables form a singlet, the 
remaining 5 = ^ 's could form either a triplet or a singlet). Thus, Haklt 
acting on \'4>vbs) gives zero. Since Haklt is a sum over projection operators, 
the lowest possible eigenvalue is zero. Hence, IV'yBs) is the ground state of 
Haklt with eigenvalue zero. The AKLT ground state (the VBS state) is spin- 
disordered and the two-spin correlation function has an exponential decay. The 
total spin of two spin-l's is 2, 1, 0. The projection operator onto spin j for a 
pair of n.n. spins has the general form 



[l{l 



t' 



^^■(^« + ^^+1) = n uLr,_.,Ju^ (43) 



where S^ = ^^ + S^^+i . For the AKLT model, j = 2 and ^ = 1, . From (42) 
and (43), 



Haklt 



E 



-{S iTS i+i) + -{S i^S i+i) +- 



(44) 



The method of construction of the AKLT Hamiltonian can be extended to higher 
spins and to dimensions d > 1. The MG Hamiltonian (apart from a numerical 
prefactor and a constant term) can be written as 
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H = Y,Pi(^^+^^+l+^^+2) (45) 

i 

The S — 1 HAFM and the AKLT chains are in the same Haldane phase, char- 
acterised by a gap in the excitation spectrum. The physical picture provided 
by the VBS ground state of the AKLT Hamiltonian holds true for real systems 
[Q . The excitation spectrum of Haklt cannot be determined exactly. Arovas 
et al. 1^^ have proposed a trial wave function 

N 

\k) = N-i J2 e'''^r \'l'yBs) ,fi = z,+,- (46) 

and obtained 

.,. {k\HvBs\k) 25 + 15co.(fc) 
'^''^ ^ {k\k) = 27 (^^^ 

The gap in the excitation spectrum A = iy at k = tt. Another equivalent way 
of creating excitations is to replace a singlet by a triplet spin configuration |^^ . 

3 Spin Ladders 

A. Undoped ladders 

In the last Section, we discussed some exact results for interacting spin sys- 
tems. The powerful technique of BA was described. The BA cannot provide 
knowledge of correlation functions. There is another powerful technique for Id 
many body systems known as bosonization ||4^ which enables one to calcu- 
late various correlation functions for Id systems. After the discovery of high- 
Tc cuprate systems, the study of 2d AFMs acquired considerable importance. 
There are, however, not many rigorous results available for 2d spin systems. 
Ladder systems interpolate between a single chain (Id) and the square lattice 
(2d) and are ideally suited for the study of the crossover from Id to 2d. Consider 
a two-chain spin ladder (Figure 4) described by the AFM Heisenberg exchange 
interaction Hamiltonian 

The n.n. intra-chain and the rung exchange interactions are of strength J and 
Jb. respectively. When Jr = 0, one obtains two decoupled AFM spin chains for 
which the excitation spectrum is known to be gapless. Dagotto et al. |Q derived 
the interesting result that the lowest excitation spectrum is separated by an 
energy gap from the ground state. The result is easy to understand in the simple 
limit in which the exchange coupling Jr along the rungs is much stronger than 
the exchange coupling J along the chains. The intra-chain coupling may thus be 
treated as perturbation. When J = , the exact ground state consists of singlets 
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along the rungs, each singlet having the spin configuration —7= [|t|,) — lit)]- The 

ground state energy is —^^^f^, where N is the number of rungs in the ladder. 
In first order perturbation theory, the correction to the ground state energy is 
zero. The ground state has total spin S = 0. A S = 1 excitation may be created 
by promoting one of the rung singlets to a S* = 1 triplet. A triplet has the spin 
configuration |TT) {S^ = +1), 75 ITi + it) (^^ = ) and Ui) {S^ = -1). A 
triplet costs an exchange energy equal to Jr . The weak coupling along the 
chains gives rise to a band of propagating S — 1 magnons with the dispersion 
relation 

uj{k) = Jr + Jcosk (49) 

in first order perturbation theory (fc is the wave vector). The spin gap, defined 
as the minimum excitation energy is given by 

AsG =Lo{ii)~ {Jr - J) (50) 

The two-spin correlations decay exponentially along the chains showing that the 
ground state is a quantum spin liquid (QSL). As the rung exchange coupling Jr 
decreases, one expects that the spin gap will also decrease and ultimately be- 
come zero at a critical value of Jr. Barnes et al. |51|, however, put forward the 
conjecture that Asg > for all ^ > 0, including the isotropic limit Jr — J. A 
variety of numerical techniques like exact diagonalization of finite-sized ladders 
[ pO| , Quantum Monte Carlo (QMC) simulations l^^l f-nd density-matrix renor- 
malization group (DMRG) ||5^ have verified the conjecture. We now consider 
the case of an n-chain ladder. A surprising fact emerging out of several theoret- 



ical studies |l54| |5^ is: the excitation spectrum has spin gap (is gapless) when n 
is even (odd). In the first (second) case, the two-spin correlation function has 
an exponential (power-law) decay. For odd n, the ladder has properties similar 
to those of a single chain. The strong coupling limit {Jr ^ J) again provides a 
physical picture as to why this is true. When n is even, the S = ^ spins along 
a rung continue to form a singlet ground state. Hence the creation of a S =1 
excitation requires a finite amount of energy as in the case of the two-chain 
ladder. The gap should decrease as n increases so that the gapless square lat- 
tice limit is reached for large n. When n is odd, each rung consists of an odd 
number of spins, each of magnitude ^ . The inter- rung (intra-chain) coupling J 
generates an effective interaction between the S = ^ rung states, which because 
of rotational invariance, should be of the Heisenberg form with an effective cou- 
pHng Jeff setting the energy scale. The equivalence of an odd-chain ladder to 
the single Heisenberg chain leads to a gapless excitation spectrum. Rojo |Bq| 
has given a rigorous proof of the gaplessness of the excitation spectrum when 
n is odd. Khveshchenko |5^ has shown that for odd-chain ladders, a topolog- 
ical term governing the dynamics at long wavelengths appears in the effective 
action, whereas, it exactly cancels for even-chain ladders. The topological term 
has similarity to the one that causes the difference between integer and half- 
odd integer spin chains. In the first case, the spin excitation spectrum has the 
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well-known Haldane gap. In the latter case, the LSM theorem shows that the 
excitation spectrum is gapless. Ghosh and Bose |Q have constructed an n-chain 
spin ladder model for which the exact ground state can be determined for all 
values of n. For n even (odd), the excitation spectrum has a gap (is gapless). 
This is true even for large n, thus the square lattice Hmit cannot be reached in 
the model. Thermodynamic properties of the S = ^ two-chain ladder have been 
first studied by Troyer et al. p9|. Using a quantum transfer matrix method, 
they obtained reliable results down to temperature T ~ 0.2J. The AFM cor- 
relation length £,AFM has been found to be 3-4 lattice spacings. The magnetic 
susceptibility x(T) shows a crossover from a Curie-weiss form, x{T) — j^at 

high temperature to an exponential fall-off, x(T) ^ '^ ^ as T ^ 0. The 
fall-off is a signature of a finite spin gap ^sG- Frischmuth et al. |Q, using a 
powerful loop algorithm, have calculated the magnetic susceptibility and found 
evidence for the gapped (gapless) excitation spectrum in the case of an even 
(odd)-chain ladder. 

A major interest in the study of ladder systems arises from the fact that 
there is a large number of experimental realizations of ladder systems. A com- 
prehensive review of major experimental systems is that by Dagotto |6J. We 
discuss here only a few interesting ladder systems. Hiroi et al. |^ were the 
first to synthesize the family of layer compounds S'r„_iCu„+i02n- Rice et al. 
[ p4| subsequently recognized that these compounds contained weakly-coupled 
ladders of ^^^^ chains. For n =3 and 5, respectively, one gets the two-chain and 
three-chain ladder compounds. Azuma et al. ||6^ have determined the tem- 
perature dependence of the magnetic susceptibility in these ladder compounds 
experimentally. A spin gap is indicated by the sharp fall of x(T)for T < SOOiiT 
in the two-chain ladder compound SrCu20^ . The magnitude of the spin gap is 
AsG ^ 420-?C. This is approximately in agreement with the theoretical result of 
A^G ~ ■^, if an exchange coupling J '-^ 1200^ is assumed. For the three-chain 
ladder compound Sr2Cu^Oz , Azuma et al. found that x(T) approaches a con- 
stant as T ^ 0, , as expected for the Id Heisenberg AFM chain. Muon spin 
relaxation measurements by Kojima et al. iQ shows the existence of a long 
range ordered state with Neel temperature Tjq = 52 K, brought about by the 
interlayer coupling. No sign of long range ordering was found in the two-chain 
ladder compound, confirming the difference between odd and even chain ladders. 
The compound LaCu02.5 is formed by an array of weakly interacting two-chain 
ladders |g^. The evidence of spin-liquid formation at intermediate tempera- 
tures (confirmed by the existence of a spin gap) and an ordered Neel state at 
low temperatures, shows that the spin singlet state is in close competition with 
a Neel state. Spin ladders, belonging to the organic family of materials, have 
also been synthesized. A recent example is the compound {C^Hi2N)2CuBr4^ 
[|6J. This compound is a good example of a strongly coupled (^ ~ 3.5 ) ladder 
system. The phase diagram of the AFM spin ladder in the presence of an ex- 
ternal magnetic field is particularly interesting. In the absence of the magnetic 
field and at T = 0, the ground state is a quantum spin liquid with a gap in the 
excitation spectrum. At a field Hc^ , there is a transition to a gapless Luttinger 
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liquid phase {g^sHc^ = /^sG , the spin gap, /is is the Bohr magneton and g the 
Lande sphtting factor) . There is another transition at an upper critical field J?C2 
to a fully polarised FM state. Both 11^ and Hc2 are quantum critical points. 
The phase transitions that occur at these points are quantum phase transitions 
as they occur at T = 0. At a quantum critical point, the system switches from 
one ground state to another. The transition is brought about by changing a 
parameter (magnetic field in the present example) other than temperature. At 
small temperatures, the behaviour of the system is determined by the crossover 
between two types of critical behaviour: quantum critical behaviour at T = 
and classical critical behaviour at T y^ 0. Quantum effects are persistent in 
the crossover region at small finite temperature and such effects can be probed 
experimentally. Refs. |^, |6^, |6^ give extensive reviews of quantum critical 
phenomena. In the case of the ladder system {CzHi2N)2GuBri^ the magneti- 
zation data, obtained experimentally, exhibit universal scaHng behaviour in the 
vicinity of the critical fields, -ffci and Hc^ ■ We remember that in the vicinity of 
critical points, physical quantities of a system exhibit scaling behaviour. Quan- 
tum spin systems provide several examples of quantum phase transitions and 
organic spin ladders are systems which provide experimental testing grounds of 
theories of such transitions. For inorganic spin ladder systems, the value of H^ 
is too high to be experimentally accessible. 
B. Frustrated spin ladders 

Bose and Gayen [Q have studied a two-chain spin ladder model with frus- 
trated diagonal couplings (Figure 5, frustrated spin systems are defined in Sec- 
tion 4) . The intra-chain and diagonal spin-spin interactions are of equal strength 
J. The exchange interactions along the rungs are of strength Jr ■ It is easy 
to show that for Jr > 2 J, the exact ground state consists of singlets along the 
rungs with the energy Eg — —^^ — where N is the number of rungs. Xian |£l| 
pointed out that the Hamiltonian of the frustrated ladder model can be written 
as 

H=JRJ2^U■'^2^ + JJ2^^■^^+^ (51) 

■i i 

where. Pi — Su + S2i , represents a composite operator at the i-th rung 
and '1' and '2' refer to the lower and upper chains respectively. Due to the 
commutativity of the rung interaction part of the Hamiltonian with the second 
part, the eigenstates of the Hamiltonian can be described in terms of the total 
spins of individual rungs. The energy eigenvalue for the state with singlets on all 
the rungs is Eg = —^^ — . The second term in the Hamiltonian (Eq.(51)) does 
not contribute in this case. If the two rung spins form a triplet, the second term 
is equivalent to the Hamiltonian of a spin-1 Heisenberg chain with a one-to-one 
correspondence between a rung of the ladder and a site of the S = 1 chain. 
Because the two parts of H commute, the eigenvalue, when the rung spins form 
a triplet is 



Ej^iJeo+'-^)N (52) 



16 



where cq = —1.40148403897(4) is the ground state energy/site of the spin 1 
Heisenberg chain. Comparing the energy Ej with the energy Eg of the rung 
singlet state, one finds that as long as ^ > (^)c = cq , the latter state is the 
exact ground state. At the critical value (^)c , there is a first order transition 
from the rung singlet state to the Haldane phase of the S = 1 chain. The 
lowest spin excitation in the rung singlet state can be created by replacing a 
rung singlet {S = 0) by a triplet {S = 1). The triplet excitation spectrum 
has no dynamics. In a more general parameter regime, i.e., when the intra- 
chain exchange interaction is not equal in strength to the diagonal exchange 
interaction, the ground and the excited states can no longer be determined 
exactly. In this case, one takes recourse to approximate analytical and numerical 
methods. Kolezhuk and Mikeska |Q have constructed a class of generaHsed 
S = ^ two-chain ladder models for which the ground state can be determined 
exactly. The Hamiltonian H is a sum over plaquette Hamiltonians and each 
plaquette Hamiltonian contains various two-spin as well as four-spin interaction 
terms. They have further introduced a toy model, the Generalised Bose-Gayen 
(GBG) model which has a rich phase diagram in which the phase boundaries 
can be determined exactly. Recently, some integrable spin ladder models with 
tunable interaction parameters have been introduced ||7^, [TJ, [T^. The integrable 
models, in general, contain multi-spin interaction terms besides two-spin terms. 
C. Doped spin ladders 

A major reason for the strong research interest in ladders is that doped 
ladder models are toy models of strongly correlated systems. The most well- 
known examples of the latter are the high-Tc cuprate systems. As already 
mentioned in the Introduction, these systems exhibit a rich phase diagram as 
a function of the dopant concentration. Doping effectively replaces the spin-^ 
's associated with the Cu^+ ions in the Cu02 planes by holes. The holes are 
mobile in a background of antiferromagnetically interacting Cu spins. Also, 
due to strong Coulomb correlations, the double occupancy of a site by two 
electrons, one with spin up and the other with spin down, is prohibited. This 
is a non-trivial many body problem because it involves a competition between 
two processes: hole delocalization and exchange energy minimization. A hole 
moving in an antiferromagnetically ordered spin background, say, the Neel state, 
gives rise to parallel spin pairs which raise the exchange interaction energy of the 
system. The questions of interest are: whether a coherent motion of the holes 
is possible, whether two holes can form a bound state (in the superconducting 
(SC) phase of the doped cuprates, charge transport occurs through the motion 
of bound pairs of holes), the development of SC correlations, the possibility 
of phase separation of holes etc. For the cuprates, a full understanding of 
many of these issues is as yet lacking (see [^ for a recent review of high-Tc 
superconductivity). The doped ladders are simple model systems in which the 
consequences of strong correlation can be studied with greater rigour than in 
the case of the structurally more complex cuprate systems. Recent experimental 
evidence ||6^ suggests that some phenomena are common to ladder and cuprate 
systems. The study of ladder systems is expected to provide insight on the 
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common origin of these phenomena. Some ladder compounds can be doped 
with holes. Much excitement was created in 1996 when the ladder compound 
Sri4-xCaxCu2404i was found to become SC under pressure at x = 13.6. The 
transition temperature Tc ~ 12K at a pressure of 3 GPa. As in the case of 
SC cuprate systems, holes form bound pairs in the SC phase of ladder systems. 
The possibility of binding of hole pairs in a two-chain ladder system was first 
pointed out by Dagotto et al. pOl. The strongly correlated doped ladder system 
is described by the t-J Hamiltonian 

Ht-j = -J2 iv (C+ Q. + H.C.) + ^ J,, 0a, - \n,n,) (53) 

(y):Ci- (ij) 

The C^and Ci^ are the electron creation and annihilation operators which act 
in the reduced Hilbert space (no double occupancy of sites). 



Cl = C+(l-n,_.) 

C^a = a,(l-n,_,) (54) 

a is the spin index and n^ , rij are the occupation numbers of the i-th and j-th 
sites respectively. The term proportional to riirij is often dropped. The first 
term in Eq.(53) describes the motion of holes with hopping integrals Ir and t for 
motion along the rung and chain respectively. In the conventional t— J ladder 
model, i and j are n.n. sites. The second term (minus the —jUirij term) is 
the usual AFM Heisenberg exchange interaction Hamiltonian. The t— J model 
thus describes the motion of holes in a background of antiferromagnetically 
interacting spins. In the undoped limit, each site of the ladder is occupied by a 
spin-i and the t— J Hamiltonian reduces to the AFM Heisenberg Hamiltonian. 
Removal of a spin creates an empty site, i.e., a hole. A large number of studies 
have been carried out on t ~ J ladder models. These are reviewed in Refs. 
[ pH 1?^ , |7q| . We describe briefiy some of the major results. A hole-doped 
single AFM chain is an example of a Luttinger Liquid (LL) which is different 
from a Fermi liquid. The latter describes interacting electron systems in higher 
dimensions and at low temperatures. A novel characteristic of a LL is spin- 
charge separation due to which the charge and spin parts of an electron (or 
hole) move with different velocities and thus become separated in space. The 
undoped two-chain ladder has a spin gap. This gap remains finite but changes 
discontinuously on doping. This is because there are now two disctinct triplet 
excitations (remember that the spin gap is the difference in energies of the lowest 
triplet excitation and the ground state). One triplet excitation is obtained 
by exciting a rung singlet to a rung triplet as in the undoped case. A new 
type of triplet excitation is obtained in the presence of two holes. A clear 
physical picture is obtained in the limit Jr ^ J. In this case, the ground state 
predominantly consists of singlets along the rungs. On the introduction of a 
hole, a singlet spin pair is broken and the hole exists with a free spin-^ . In the 
presence of two holes on two separate rungs, the two free spin-i 's combine to 
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give rise to an excited triplet state. The ground state is a singlet and consists 
of a bound pair of holes. The binding of holes can be understood in a simple 
manner. Two holes located on two different rungs break two rung singlets and 
the exchange interaction energy associated with two rungs is lost in the process. 
If the holes are located on the same rung, the exchange interaction energy of only 
one rung is lost. If Jr is much greater than the other parameters of the system, 
the holes preferentially occupy rungs in pairs. As Jr decreases in strength, the 
hole bound pair has a greater spatial extent. The lightly doped ladder system 
is not in the LL phase, i.e., no spin-charge separation occurs. The system is in 
the so-called Luther-Emery phase with gapless charge excitations and gapped 
spin excitations. A variety of numerical studies show that the hole pairs and the 
spin gap are present even in the isotropic limit Jr — J . Also, the relative state 
of hole pairs has approximate "d-wave" symmetry with the pairing amplitude 
having opposite signs along the rungs and the chains. The d-wave symmetry 
is a feature of strong correlation and is considered to be the symmetry of the 
pairing state in the case of cuprate systems. 

Bose and Gayen ||70| , [79| , pO| have constructed a two-chain t— J ladder model 
with frustrated diagonal couplings. The intra-chain n.n. and the diagonal hop- 
ping integrals have the same strength t. The other parameters have been defined 
earher. The special structure of the model enables one to determine the exact 
ground and excited states in the cases of one and two holes. The most signif- 
icant result is an exact, analytic solution of the eigenvalue problem associated 
with two holes in the infinite t — J ladder. The binding of holes has been ex- 
pHcitly demonstrated and the existence of the Luther-Emery phase established. 
For conventional t — J ladders (the diagonal bonds are missing), the only ex- 
act results that have been obtained are through numerical diagonalization of 
finite-sized ladders. Derivation of exact, analytical results in this case has not 
been possible so far. The reason for this is that as a hole moves in the anti- 
ferromagnetically interacting spin background, spin excitations in the form of 
parallel spin pairs are generated. Proliferation of states with spin excitations 
makes the solution of the eigenvalue problem extremely difficult. In the case of 
the frustrated t — J ladder model, there is an exact cancellation of the terms 
containing parallel spin pairs ^^. Thus the hole has a perfect coherent motion 
through the spin background. Frahm and Kundu |gl| has constructed an inte- 
grable t— J ladder model and obtained the phase diagram. The model contains 
terms describing correlated hole hopping in chains which may not be realizable 
in real systems. Several studies have been carried out on the two-chain Hubbard 
ladder as well as on multi-chain Hubbard and t— J ladders. References of some 
of the studies may be obtained from ||6l|| . 

4 Frustrated spin models in 2d 

In Sections 2 and 3 we have discussed quasi-ld interacting spin systems, namely, 
spin chains and ladders. As already mentioned in the Introduction, the Cu02 
plane of the undoped cuprate systems is a 2d AFM. The undoped cuprates 
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exhibit AFM LRO below a Neel temperature Tjv . On the introduction of a 
few percent of holes, the AFM LRO is rapidly destroyed leaving behind spin- 
disordered states in the CuO^ planes. This fact has triggered lots of interest in 
the study of spin systems with spin disordered states as ground states. Frus- 
trated spin models are ideal candidates for such systems. To understand the 
origin of frustration, consider the AFM Ising model on the triangular lattice. 
An elementary plaquette of the lattice is a triangle. The Ising spin variables 
have two possible values, ±1, corresponding to up and down spin orientations. 
An antiparallel spin pair has the lowest interaction energy — J . A parallel spin 
pair has the energy -I- J . In an elementary triangular plaquette, there are three 
interacting spin pairs. Due to the topology of the plaquette, all the three pairs 
cannot be simultaneously antiparallel. There is bound to be at least one paral- 
lel spin pair. The parallel spin pair may be located along any one of the three 
bonds in the plaquette and so the ground state is triply degenerate. The Ising 
model on the full triangular lattice has a highly degenerate ground state such 
that the entropy/spin is a finite quantity. As a result, the system never orders 
including at T = . Frustration occurs in the system since all the spin pair 
interaction energies cannot be simultaneously minimised. On the other hand, 
consider the AFM Ising model on the square lattice. All the four spin pairs 
in an elementary square plaquette can be made antiparallel and so there is no 
frustration. The system exhibits magnetic order below a critical temperature. 
If one of the spin pair interactions in each elementary square plaquette is FM 
and the rest AFM, frustration occurs in the square lattice spin system. A spin 
system with mixed FM and AFM interactions is frustrated if the sign of the 
product of exchange interactions around an elementary plaquette is negative. 
In the case of a purely AFM model, frustration occurs if the number of bonds 
in an elementary plaquette of the lattice is odd. Examples of such lattices in 2d 
are the triangular and kagome lattices. In 3d, the pyrochlore lattice, the elemen- 
tary plaquette of which is a tetrahedron provides an example. A spin system 
is also frustrated due to the presence of both n.n. as well as further-neighbour 
interactions. Consider AFM n.n. as well as n.n.n. interactions between a row 
of three Ising spins. Again, all the three spin pairs cannot simultaneously be 
made antiparallel. 

Let us now treat the spins as classical vectors {S -^ oo). For AFM spin-spin 
interaction, the lowest energy is achieved for an antiparallel spin configuration. 
In the classical limit, the spins on a bipartite lattice are ordered in the AFM 
Neel state. On a non-bipartite lattice, such as the triangular lattice, the classical 
ground state represents a compromise between competing requirements. In the 
ground state, the spins form an ordered three-sublattice structure with 120" 
between n.n. spins on different sublattices. The ground state of the classical 
Heisenberg model on the kagome lattice is, however, highly degenerate and 
disordered. We now consider the full quantum mechanical spin Hamiltonian and 
ask the question how the classical ground states are modified when quantum 
fiuctuations are taken into account. In the case of the triangular lattice, it is now 
believed that the quantum mechanical ground state of the S — ^ HAFM model 
has AFM LRO of the Neel-type, i.e., quantum fiuctuations do not destroy the 
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three-sublattice order of the classical ground state. In the second scenario, when 
the classical ground state is highly degenerate and disordered, thermal/quantum 
fluctuations select a subset of states which tend to incorporate some degree of 
long range order. This is the phenomenon of 'order from disorder' which is 
counterintuitive since order is brought about by fluctuations which normally 
have disordering effects. The classical kagome lattice HAFM ground states 
include both coplanar as well as noncoplanar spin arrangements and fluctuations 
lead to the selection of coplanar order. This kind of ordering is particularly true 
for large values of the spin S. As the magnitude of the spin is decreased towards 
S = ^, the quantum fluctuations increase in strength. These fluctuations often 
destroy the ordered structure obtained for large S. The quantum mechanical 
ground states of the S = ^ HAFM on the kagome and pyrochlore lattices have 
been found to be spin disordered. Some recent references of frustrated magnetic 
systems are |g^, |8^ . The triangular lattice S = ^ HAFM is the flrst example 
of a spin model in which frustration occurs due to lattice topology JSJ]. The 
5 = i HAFM model has also been studied on a partially frustrated pentagonal 
lattice [?] and a parameter region identifled in which the ground state has AFM 
LRO of the Neel-type. 

Two well-known examples of spin-disordered states are the quantum spin 
Hquid (QSL) and dimer or valence bond (VB) states. A QSL state is a spin 
singlet with total spin S = and has both spin rotational and translational 
symmetry. In a VB state, spin rotational symmetry is present but tanslational 
symmetry is broken. In such states pairs of spins form singlets which are called 
VBs or dimers with the VBs being frozen in space. A well-known example of a 
QSL state is the resonating-valence-bond (RVB) state [Q which is a coherent 
Hnear superposition of VB states (Figure 6). The RVB state is the starting 
point for the well-known RVB theory of high- Tc SC. Spin-disordered (no AFM 
LRO as defined in Eq. (8)) states with novel order parameters are: 

(a) Chiral states 

In these states, the spins are arranged in configurations characterised by the 
order parameter 

^^ = {'^U'^.+7^^.+y) (55) 

with the three spins belonging to one plaquette of the square lattice and x, y 
denoting unit vectors in the x and y directions respectively. The chiral state 
breaks time reversal symmetry or a refiection about an axis (parity). 

(b) Dimer states 

These are the VB states in which the VBs are frozen in space. A well-known 
example of such states is the columnar dimer (CD) states. In such states, 
the VBs are arranged in columns. On the square lattice, four such states are 
possible. The order parameter of CD states is 

A = Ul)^l.0,^- + ^^,,_-~^,--^^,-)) (56) 



l-\-x l+y l — x ^ — y^ 

where the 1-sites are even and rj{l) = +1(— 1) if both Ix and ly are even (odd). 
The order parameter takes the values 1, i, —1, —i for the four CD states shown 
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in Figure 3. 

(c) Twisted states: 

At the classical level {S — > cxd), the spins in the twisted state are arranged 
in incommensurate structures. These configurations can be visualised as spins 
lying in a plane and with a twist angle in some direction. It is possible that 
such states survive the inclusion of quantum fluctuations. The order parameter 
is vectorial in nature and is given by 

Ti = {~^i^C^,+7 + ^i+y)) (57) 

These states are called spin nematics and are different from helimagnets in which 
both Ti and the spin-spin correlation functions show LRO. 

(d) Strip or Collinear States 

In a classical picture, the spins are ferromagnetically ordered in the x direction 
and antiferromagnetically ordered in the y direction. The configuration obtained 
by rotating the previous one by ^ is also possible. The order parameter is given 
by 

Ci = ( 5^^(5^,+^- ^i+y + ^i-^-^i-y)) (58) 

Ci takes the values 1, —1 in the two different strip configurations. 

The spin-disordered states decribed above are quantum-coherent states and 
are characterised by novel order parameters. The term 'quantum paramagnet' 
is often used to describe such states. 

Examples of real frustrated systems are many [g^i Hi ■ The best studied ex- 
perimental kagome system is the magnetoplumbite, SrCrs-xGa4+xOig . The 
system consists of dense kagome layers of 5 = | Cr ions, separated by dilute tri- 
angular layers of Cr. In a mean-field theory of the HAFM, the high temperature 
susceptibility is given by 

-t + Ocw 

2 2 

Here, the Curie constant C = ^^^ , where ^b is the Bohr magneton and, 
p — g[s{s -t- 1)]2, g being the Lande splitting factor governing the splitting of 
the spin multiplet in a magnetic field. Also, 6cw is the Curie- Weiss temperature. 
The Neel ordering temperature T/v is defined experimentally using bulk probes. 
One looks for singularities in either the specific heat C(T) or the temperature 
derivative of the susceptibility x(r) . In the case of non-frustrated systems, 
Tn ^ dew and in the second case, T/v << Ocw ■ Since a frustrated system 
may not order at all in a conventional sense, the hallmark of such a system is 
Tc << 9cw where Tc is the temperature below which new types of spin order set 
in. The Curie-weiss temperature 9cw is an experimentally measurable quantity. 
One defines an empirical measure of frustration by the quantity 

/ = -|^ (60) 
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Frustration corresponds to / > 1 . For the kagome AFM SrCrsGa^Oig , f is as 
high as 150. 

SrCrGaO displays unconventional low-T behaviour ||8^. One of these is 
the insensitivity of the specific heat C{T) to applied magnetic fields H as large 
as twice the temperature. Susceptibility measurements show the existence of a 

gap AsG in the triplet spin excitation spectrum. For T < Asg,X ^ e ''bt . 
The specific heat, however, does not decrease exponentially for T < Asg, i.e., 
does not have a thermally activated behaviour. It has a T^ dependence. This 
fact along with the experimental observation of insensitivity of C{T) to external 
magnetic field have been explained by suggesting that a large number of singlet 
excitations fall within the triplet gap |^^ . Numerical evidence of such excita- 
tions has been obtained in the case of the S = ^ HAFM on the kagome lattice 
[ p7| . The number of such excitations has been found to be ~ (1.15)^ where 
N is the number of spins in the lattice. Mambrini and Mila [SSJ have recently 
established that a subset of short-range RVB states captures the specific low 
energy physics of the kagome lattice HAFM and the number of singlet states 
in the singlet-triplet gap is (1.15)^ in agreement with the numerical results. 
The appearance of singlet states in the singlet-triplet gap could be a generic 
feature of strongly frustrated magnets. Other examples of such systems are: 
the S = ^ frustrated HAFM on the ^-depleted square lattice describing the 2d 



AFM compound CaV^Og ||89|, gO|, the HAFM on the 3d pyrochlore lattice and 



a Id system of coupled tetrahedra [^ . Bose and Ghosh ||90|| have constructed a 
frustrated S = ^ AFM model on the ^—depleted square lattice and have shown 
that in different parameter regimes the plaquette RVB (PRVB) and the dimer 
states are the exact ground states. In the PRVB state, the four-spin plaquettes 
(Figure 7) are in a RVB spin configuration which is a linear superposition of 
two VB states. In one such state, the VBs (spin singlets) are horizontal and 
in the other state the VBs are vertical. In the dimer state, VBs or dimers 
form along the bonds joining the four-spin plaquettes. Both the PRVB and 
the dimer states are spin disordered states. The state intermediate between the 
PRVB and the dimer states has AFM LRO. Both the PRVB and dimer phases 
are characterised by spin gaps in the excitation spectrum. For the unfrustrated 
HAFM model on the ^ -depleted square lattice, Troyer et al. ||9^ have carried 
out a detailed study of the quantum phase transition from an ordered to a dis- 
ordered phase. In the ordered phase, the excitation spectrum is gapless. The 
spin gap AsG continuously goes to zero in a power-law fashion at the quan- 
tum critical point separating a gapped disordered phase from a gapless ordered 
phase. Chung et al. |^^ have carried out an extensive study on the possible 
paramagnetic phases of the Shastry-Sutherland model [^. In addition to the 
usual dimer phase, they find the existence of a phase with plaquette order and 
also a topologically ordered phase with deconfined S = ^ spinous and heHcal 
spin correlations. Takushima et al. ||9j| have studied the ground state phase 
diagram of a frustrated S = ^ quantum spin model on the square lattice. This 
model includes both the Shastry-Sutherland model as well as the spin model on 
the ^—depleted square lattice as special cases. The nature of quantum phase 
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transitions among the various spin gap phases and the magnetically ordered 
phases has been clarified. 

Lieb and Schupp [Q have derived some exact results for the fully frustrated 
HAFM on a pyrochlore checkerboard lattice in 2d. This lattice is a 2d version 
of the 3d pyrochlore lattice. The ground states have been rigorously proved to 
be singlets. The Lieb-Mattis theorem is applicable only for bipartite lattices. 
Hence, the proof for the pyrochlore checkerboard lattice is a new result. Lieb and 
Schupp have further proved that the magnetization in zero external field vanishes 
separately for each frustrated tetrahedral unit. Also, the upper bound on the 
susceptibility is | in natural units for both T — and T 7^ . Frustration can 
also occur from a competition between exchange anisotropy and the transverse 
field terms as in the case of the transverse Ising model. Moessner et al. |^g| 
have shown that the transverse Ising model on the triangular (kagome) lattice 
has an ordered (disordered) ground state. 

Frustrated AFMs with short-range dimer or RVB states as ground states 
have a gap in the spin excitation spectrum and the two-spin correlation func- 
tion has an exponential decay as a function of the distance separating the spins. 
A single branch describes the S' = 1 triplet excitation spectrum. In Section 2 
we pointed out that in the case of the S = ^ HAFM chain, a pair of spinous 
(each spinon has spin S = ^ ) are the fundamental excitations. The lowest 
excitation spectrum is thus not a single branch of 5 = 1 magnon excitations 
but a continuum of scattering states with well-defined lower and upper bound- 
aries. AFM compounds in 2d, in general, have excitation spectra described 
by S* = 1 magnons. Anderson |Q suggested that a RVB state may support 
pairs of spinous as excitations which are deconfined via a rearrangement of the 
VBs. In this case, an extended and highly dispersive continuum of excitations 
is expected. Recently, Coldea et al. |97J have investigated the ground state or- 
dering and dynamics of the 2d S" = ^ frustrated AFM CS2CUCI4 using neutron 
scattering in high magnetic fields. The dynamic correlations exhibit a highly 
dispersive continuum of excited states which are characteristic of the RVB state 
and arise from pairs of S" = ^ spinous. A recent paper 'RVB Revisited' by 
P.W.Anderson [Q focuses on the relevance of the RVB state to describe the 
normal state of the Cu02 planes in the high- Tc cuprate superconductors. 

5 Concluding remarks 

In Sections 1-4, a brief overview of low-dimensional quantum magnets, specially, 
antiferromagnets has been given. The subject of quantum magnetism has wit- 
nessed an unprecedented growth in research activity in the last decade. This is 
one of the few research areas in which rigorous theories can be worked out and 
experimental realizations are not diffficult to find. Coordination chemists and 
material scientists have prepared novel materials and constructed new applica- 
tion devices. Experimentalists have employed experimental probes of all kinds 
to refine old data and uncover new phenomena. Theorists have taken recourse 
to a variety of analytical and numerical techniques to explain the experimen- 
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tally observed properties as well as to make new predictions. Powerful theorems 
have been proposed and exact results obtained. This trend is still continuing 
and will lead to further breakthroughs in materials, phenomena and techniques 
in the coming years. The study of doped magnetic materials which include 
cuprates, ladders and spin chains has acquired considerable importance in re- 
cent times. The dopants can be magnetic and nonmagnetic impurities as well as 
holes. To give one example, it has been possible to dope the spin-1 Haldane-gap 
AFM compound Y2BaNiO^ with holes. Neutron scattering experiments reveal 
the existence of midgap states and an incommensurate double-peaked structure 
factor [Q. Several new experiment al re sults on the Haldane-gap AFMs have 



been obtained in the last few years | 100 | which add to the richness of phenom^ 



ena observed in magnetic syatems. In this overview, only a few of the aspects 
of quantum magnetism have been highlighted. Conventional 2d and 3d mag- 
nets have not been discussed at all as a good understanding of these materials 
already exists. We have not discussed recent advances in material appHcations 
which include materials exhibiting colossal magnetoresistance, molecular mag- 
nets, nanocrystalline magnetic materials, magnetoelectronic devices (the study 
of which constitutes the new subject of spintronics) etc. A report on some of 
these developments as well as some recent issues in quantum magnetism may 



be obtained from Ref. [|101 | 
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Figure Captions 

Figure 1. The Shastry-Sutherland model 

Figure 2. Five types of interaction in the Ji — J2 — J3 — J4 — J5 model 

Figure 3. Four columnar dimer states 

Figure 4. A two-chain spin ladder 

Figure 5. The two-chain frustrated spin ladder model 

Figure 6. An example of a RVB state 

Figure 7. The 4— depleted square lattice. 
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